Unifying inflation with dark energy in modified F(R) Hofava-Lifshitz gravity 
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We study FRW cosmology for a non-linear modified F(R) Hofava-Lifshitz gravity which has a 
viable convenient counterpart. A unified description of early-time inflation and late-time acceleration 
is possible in this theory, but the cosmological dynamic details are generically different from the 
ones of the convenient viable F(R) model. Remarkably, for some specific choice of parameters they 
do coincide. The emergence of finite-time future singularities is investigated in detail. It is shown 
that these singularities can be cured by adding an extra, higher-derivative term, which turns out to 
be qualitatively different when compared with the corresponding one of the convenient F(R) theory. 

PACS numbers: 98.80.-k,04.50.+h,11.10.Wx 



I. INTRODUCTION 



Current observational data clearly indicates that our universe has undergone at least two periods of accelerated 
expansion: the early-time inflation and the present late-time cosmic acceleration. In spite of the existence of a 
number of (partially) successful scenarios for the inflationary and dark energy epochs, the fundamental issue of a 
unified description of the whole cosmic history scenario remains open. One possibility to solve this problem, relying 
only on the presence of gravity (and not on the introduction of additional cosmological fields), is modified gravity 
(for a quick presentation, see [l(). Indeed, this approach suggests a very natural unification of early-time inflation 
and late-time cosmic acceleration, as a purely gravitational alternative (see Q for a review of such unified models of 
f^i ' modified gravity). But, how general is this scenario? 

The Hofava-Lifshitz quantum gravity @ has been conjectured to be renormalizable in four dimensions, at the 
price of explicitly breaking Lorentz invariance. Its generalization to an i 7 '(i?)-formulation, which seems to be also 
renormalizable in 3 + 1 dimensions, has been considered in Refs. [1, @, where the Hamiltonian structure and FRW 
cosmology, in a power-law theory, have been investigated for such modified F(R) Hofava-Lifshitz gravity. It was also 
conjectured there that it sustains, in principle, the possibility of a unified description of early-time inflation and the 
dark energy epochs. 

The purpose of the present work is to study a realistic non-linear F(R) gravity in the Hofava-Lifshitz formulation, 
with the aim to understand if such theory is in fact directly able to predict in a natural way the unification of the 
two acceleration eras, similarly as its done in the convenient version. It has been shown in Q that, for a special 
choice of parameters, the FRW equations do coincide with the ones for the related, convenient F(R) gravity. This 
means, in particular, that the cosmological history of such Hofava-Lifshitz F(R) gravity will be just the same as for 
its convenient version (whereas black hole solutions are generically speaking different). For the general version of the 
theory the situation turns out to be more complicated. Nevertheless, the unification of inflation with dark energy is 
still possible and all local tests can also be passed, as we will prove. 

Note that in the Hofava-Lifshitz formulation of F(R) gravity, the Lorentz symmetry is explicitly broken and the 
restoration of the Lorentz symmetry at the observed energy scale is the main problem in this formulation. The 
obtained metrics in this paper are, however, FRW metrics which are almost Lorentzian at the scale of galaxy or solar 
system. Then when the matter sector has a Lorentz symmetry in the flat background, if the Newton law is reproduced, 
the violation of the Lorentz symmetry could be difficult to be observed. The restoration of the Newton law is actively 
investigated but in this paper, we just show a mechanism that the scalar field corresponding to the so-called scalaron 
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in the usual F(i?)-gravity does not give an observable correction to the Newton law. Although there is a problem 
with the Lorentz symmetry breaking, the Hofava-Lifshitz F(R) gravity has a much richer structure than standard 
F(R) gravity. 

The paper is organized as follows. Next section briefly reviews the formulation of modified F(R) Hofava-Lifshitz 
gravity and its corresponding FRW cosmology. The reconstruction of the theory is presented in Sect. IILTI It is demon- 
strated that different functional forms of the theory may reproduce the same ACDM era, phantom (super) acceleration, 
or any other given cosmology. Sect. IIVI is devoted to the analysis of F(R) Hofava-Lifshitz theories whose convenient 
counterparts have been proposed as viable candidates for inflation-dark energy unification. The unification is again 
possible but with rather different physical properties. Moreover, corrections to Newton's law are negligible in the mod- 
els under consideration (Sect. [Vj) . The emergence of finite-time future singularities and their avoidance is discussed 
in Sect. I VII In particular, it turns out that even when the viable modified gravity is non-singular, its Hofava-Lifshitz 
counterpart may still remain a singular theory. Some discussions and an outlook are provided in the last section. 



II. MODIFIED F(R) HORAVA-LIFSHITZ GRAVITY 

In this section, modified Hofava-Lifshitz F(R) gravity is briefly reviewed [4, 5]. We start by writing a general metric 
in the so-called ADM decomposition in a 3 + 1 spacetime (for more details see @, Q and references therein), 

ds 2 = ~N 2 dt 2 + glf (dx l + N l dt)(dx j + N j dt) , (I) 

where i,j = 1, 2, 3, N is the so-called lapse variable, and N l is the shift 3-vector. In standard general relativity (GR), 
the Ricci scalar can be written in terms of this metric, and yields 

R = K l3 K lj - K 2 + R.^ + 2V^V v n v - n»V„n>*) , (2) 

here K = g l i Kij, K~ij is the extrinsic curvature, R^ is the spatial scalar curvature, and a unit vector perpendicular 
to a hypersurface of constant time. The extrinsic curvature K%j is defined as 

K<i = 7^7 (dif - Vf ] N, - V^nA . (3) 



LJ 3 27V V -' •? 

In the original model Q, the lapse variable N is taken to be just time-dependent, so that the projectability condition 
holds and by using the foliation- preserving diffeomorphisms ([5]), it can be fixed to be N = 1. As pointed out in 31], 
imposing the projectability condition may cause problems with Newton's law in the Hofava gravity. On the other 
hand, Hamiltonian analysis shows that the non-projectable F(R)-mode\ is inconsistent fllj . For the non-pro jectable 
case, the Newton law could be restored (while keeping stability) by the "healthy" extension of the original Hofava 
gravity of Ref . [3l| . 

The action for standard F(R) gravity can be written as 

S = J d 4 xy/g~WNF(R) . (4) 

Gravity of Ref. Q is assumed to have different scaling properties of the space and time coordinates 

x 1 -> bx 1 , t-> b z t , (5) 

where z is a dynamical critical exponent that renders the theory renormalizable for z — 3 in 3 + 1 spacetime dimensions 
Q (For a proposal of covariant renormalizable gravity with dynamical Lorentz symmetry breaking, see Q). GR is 
recovered when z = 1. The scaling properties ([5]) render the theory invariant only under the so-called foliation- 
preserving diffeomorphisms: 

8x l = Q{x\t), St = f{t). (6) 

It has been pointed that, in the IR limit, the full diffeomorphisms are recovered, although the mechanism for this 
transition is not physically clear. The action considered here was introduced in Ref. [4|, 

S = -L / dtd 3 xy^)NF(R) , R = K tJ K %1 - \K 2 + i? (3) + 2/xV 'Jn^V ' v n v - n^n") - i (3) (<?[f) » (7) 
2k z J J 

where k is the dimensionless gravitational coupling, and where, two new constants A and [i appear, which account for 
the violation of the full diffeomorphism transformations. A degenerate version of the above F"(i?)-theory with fi = 
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has been proposed and studied in Ref. @ . Note that in the original Hof ava gravity theory Q , the third term in the 
expression for R can be omitted, as it becomes a total derivative. The term L^ 3 \g^) is chosen to be Q 

LW(9% } )=EVG ijkl E kl , (8) 
where Gijki is the inverse of the generalized De Witt metric, namely 

Qijkl _ 1 ^(3) tk g (3) jl + g (3) U g (3) jk^j _ Xg (3) ij g {3) kl _ ^ 



Therefore we have 

Gijki 



2 \9ik 9ji + 9a 9jk J ~ A 9ij 9u > A - 3A _ 1 ■ l iL U 



Note that G l ^ kl is singular for A = 1/3 and therefore G^u exist if A ^ 1/3. 

In Ref. 0], the expression for Eij is constructed to satisfy the "detailed balance principle" in order to restrict the 
number of free parameters of the theory. This is defined through variation of an action 

v^' = ^#, (ID 

where the form of W[gjy^] is given in Ref. [l(| for z = 2 and z — 3. Other forms for L^(g^') have been suggested 
that abandons the detailed balance condition but still render the theory power-counting renormalizable (see Ref. 0]). 

We are interested in the study of (accelerating) cosmological solutions for the theory described by action (|7J|. 
Spatially-flat FRW metric is assumed 



ds 2 = ~N 2 dt 2 + a 2 {t) { dxl f ■ (12) 
»=i 



If we also assume the projectability condition, N can be taken to be just time-dependent and, by using the foliation- 
preserving diffeomorphisms ([5]), it can be fixed to be unity, N = 1. When we do not assume the projectability 
condition, N depends on both the time and spatial coordinates, first. Then, just as an assumption of the solution, N 
is taken to be unity. 

For the metric (fT2j) . the scalar R is given by 

3(l-3A + 6^)ff 2 6^ d ( H 



R = ^ W^ + fdt{N)- (13) 

For the action (J7J, and assuming the FRW metric (fTB"]) . the second FRW equation can be obtained by varying the 

(3) 

action with respect to the spatial metric g^- , which yields 

= F(R) - 2(1 - 3A + 3/i) (H + 3i? 2 ) F'(R) - 2(1 - 3\)RF"{R) + 2/i (r 2 F {3) (R) + RF"{R)j + n 2 Pm , (14) 

here k 2 = 16irG, p m is the pressure of a perfect fluid that fills the Universe, and N = 1. Note that this equation 
becomes the usual second FRW equation for convenient F(R) gravity by setting the constants A = /i = 1. When 
we assume the projectability condition, variation over N of the action (JT)) yields the following global constraint 







J d 3 x F{R) - 6(1 - 3A + 3^)H 2 - 6pH + &fj,HRF"(R) - n 2 p m . (15) 



Now, using the ordinary conservation equation for the matter fluid p m + 3H(p m + p m ) = 0, and integrating Eq. (|14l) . 

= F(R) - 6 [(1 - 3A + 3(i)H 2 + ak] F'(R) + %pHRF" (R) - K 2 p m - , (16) 

L J o J 

where C is an integration constant, taken to be zero, according to the constraint equation (|15[) . If we do not assume 
the projectability condition, we can directly obtain (jTHJ) , which corresponds to the first FRW equation, by variation 
over N. Hence, starting from a given F(R) function, and solving Eqs. (fT4"|) and (fl"5|) . a cosmological solution can be 
obtained. 
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III. RECONSTRUCTING FRW COSMOLOGY IN F(R) HORAVA-LIFSHITZ GRAVITY 

To start, let us analyze the simple model F(R) = R, which cosmology was studied in [13| (for a complete analysis 
of cosmological perturbations, see [3). In such a case, the FRW equations look similar to GR, 

H2 = W^i) Pm ' * = -whT) {Pm+Pmh (17) 

where, for A — y 1, the standard FRW equations are recovered. Note that the constant p is now irrelevant because, as 
pointed out above, the term in front of p in ([7]) becomes a total derivative. For such theory, one has to introduce a 
dark energy source as well as an inflaton field, in order to reproduce the cosmic and inflationary accelerated epochs, 
respectively. It is also important to note that, for this case, the coupling constant is restricted to be A > 1/3, otherwise 
Eqs. (fl~7|) become inconsistent. It seems reasonable to think that, for the current epoch, where R has a small value, 
the IR limit of the theory is satisfied A ~ 1, but for the inflationary epoch, when the scalar curvature R goes to 
infinity, A will take a different value. It has been realized that, for A = 1/3, the theory develops an anisotropic Weyl 
invariance (see Q), and thus it takes an special role, although for the present model this value is not allowed. 

We now discuss some cosmological solutions of F(R) Hofava-Lifshitz gravity. The first FRW equation, given by 
(ITS)) with C — 0, can be rewritten as a function of the number of e-foldings 77 = In ^ , instead of the usual time t. This 
technique has been developed in Ref. [l5| for convenient F(R) gravity (for review of reconstruction method in modified 
gravity, see [l6j]), where it was shown that any F(R) theory can be reconstructed for a given cosmological solution. 
Here, we extend such formalism to the Hofava-Lifshitz F(R) gravity. Since = H-jj^ and — H 2 ' + H^-^, 
the first FRW equation (IT5|) is rewritten as 

= F(R) - 6 [(1 - 3A + 3p)H 2 + pHH 1 ] + 36pH 2 [(1 - 3A + Qp)HH' + pH' 2 + pH"H] - F } R ^ - K 2 p m , 

dR d~R 

(18) 

where the primes denote derivatives with respect to 77. Thus, in this case there is no restriction on the values of A or 
p. By using the energy conservation equation, and assuming a perfect fluid with equation of state (EoS) p m — w rn p 
the energy density yields 



m . 



Pm = P0 a-^ + ^ = p af {1+w ^e-^ 1+w ^ . (19) 

As the Hubble parameter can be written as a function of the number of e-foldings, H = H(rj), the scalar curvature in 
(|T3"| takes the form 

R = 3(1 - 3A + 6p)H 2 + 6pHH' , (20) 

which can be solved with respect to 7/ as 77 = rj(R), and one gets an expression (fl8|) that gives an equation on F(R) 
with the variable R. This can be simplified a bit by writing G(rj) = H 2 instead of the Hubble parameter. In such 
case, the differential equation (TT51) yields 

= F(R) - 6 [(1 - 3A + 3fx)G + £g'1 + 18/z [(1 - 3A + 6 M )GG' + pGG"} - n 2 Po a^ (1+w) e~^ 1+w ^ , 

L 2 J dR d?R 

(21) 

and the scalar curvature is now written as R = 3(1 — 3A + 6/i)G + 3/iG'. Hence, for a given cosmological solution 
H 2 = G(rj), one can resolve Eq. (|2"Tj) . and the F(R) that reproduces such solution is obtained. 

As an example, we consider the Hubble parameter that reproduces the ACDM epoch. It is expressed as 

H 2 = G(r7) = H 2 + ^ Po a- 3 = H 2 + yAVe' 3 " . (22) 

where Hq and po are constant. In General Relativity, the terms on the rhs of Eq. (|2"2"|) correspond to an effective 
cosmological constant A = 3Hq and to cold dark matter with EoS parameter w — 0. The corresponding F(R) 
can be reconstructed by following the same steps as described above. Using the expression for the scalar curvature 
R = 3(1 — 3A + 6/i)G + 3/iG', the relation between R and r\ is obtained, 

3J , _ £- 3 (l-3A + 6 M )ff 2 

3fc(l + 3(/x - A)) ' 1 ' 
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where k = ^pQa 3 . Then, substituting (|22|) and (|23|) into Eq. (|2~T|) , one gets the differential expression 



= (1 - 3A + 3fi)F(R) -2[l-3X+-n]R+ 9/x(l - 3A)iJ 2 



-6p(R - 9(J,H$)(R - 3H 2 (1 - 3A + 6/i)) - R - 3(1 - 3A + 6fi)H 2 , (24) 

where, for simplicity, we have considered a pressureless fluid w = in Eq. (|21[) . Performing the change of variable 
x = 3 H'^(i+3{p,-\)) ' ^ nG homogeneous part of Eq. (|24[) can be easily identified as an hypergeometric differential equation 

d 2 F dF 
= x(l - x) — + (7 - {a + j3 + l)x) — - apF , (25) 
dx A dx 

with the set of parameters (a, /?, 7) being given by 

The complete solution of Eq. ([25]) is a Gauss' hypergeometric function plus a linear term and a cosmological constant 
coming from the particular solution of Eq. (I24[) . namely 

F(R) = &F(a, /?, 7; x) + C 2 x 1 ~~<F{a - 7 + l,/3 - 7 + 1, 2 - 7; x) + —R - 2A . (27) 

where C\ and C2 are constants, K\ = 3 A — 1 and A = — 3 ^°i_3A-f^y ■ Note that for the exact cosmology (|22|) . the 
classical F(R) gravity was reconstructed and studied in Refs. [15j-[l7J. In this case, the solution (1271) behaves similarly 
to the classical F(R) theory, except that now the parameters of the theory depend on (A,/i), which are allowed to 
vary as it was noted above. We can conclude that the cosmic evolution described by the Hubble parameter (|22p is 
reproduced by this class of theories. 

One can also explore the solution (|22[) for a particular choice on the parameters fi = A — i, which plays a special 
role as it is shown below. In this case, the scalar R turns out to be a constant, and Eq. ((24j) , in the presence of a 
pressureless fluid, has the solution 

1 - 3 
F(R) = — R - 2A , with A=-(3\-l)Hg. (28) 

Kl 2 

Hence, for this constraint on the parameters, the only consistent solution reduces to the Hofava linear theory with a 
cosmological constant. 

As a further example, we consider the so-called phantom accelerating expansion. Currently, observational data do 
not totally exclude the possibility that the Universe could have already crossed the phantom divide, which means 
that the effective EoS for dark energy would presently be slightly less than — 1. Such kind of system can be easily 
expressed in GR, where the FRW equation reads H = ^Pph- Here the subscript "ph" denotes the phantom nature 
of the fluid, which has an EoS given by p p h = w p hP P h with w p h < — 1. By using the energy conservation equation, 
the solution for the Hubble parameter turns out to be 

H(t) = ^- t , (29) 

where Hq — —1/3(1 + iu p h), and t s is the Rip time which represents the time still remaining up to the Big Rip 
singularity. As in the above example, one can rewrite the Hubble parameter as a function of the number of e-foldings; 
this yields 

G(r)) = H 2 ( V ) = H 2 e 2, i/ Ho . (30) 

Then, by using the expression of the scalar curvature, the relation between R and ry is given by 

e 2n/Ho = R (31) 

H (AH + 6fj) ' V ' 
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By inserting and (j3"Tj) into the differential equation (|2"Tj) . we get 

di? 2 dR 



k2 <PF{R)_ + kiR dF(R± + koF(R) = Q (32) 



(Ag + 6 M )((Ag + 3 M )) (Aff + 6/i) 2 

1 6/x(AHq + 12/x) ' ° 12/x(Aff + 12/x) ' 



where 



here we have neglected any kind of matter contribution for simplicity. Eq. (1321) is an Euler equation, whose solution 
is well known 

F(R) = C 1 R™ + +C 2 R™-, where m ± = Llh±VEE^E^h . (34 ) 

Such theory belongs to the class of models with positive and negative powers of the curvature introduced in Ref. [l8[ . 
The existence of the negative power of curvature indicates that the flat Minkowski space is not realized. If the 
coefficient C2 is small enough, however, the difference from the flat Minkowski space could not be observed. Moreover, 
one of the reasons to consider non-linear models of next section is related with the fact that Minkowski solution may 
be realized there. 

Hence, a F(R) Hofava-Lifshitz gravity has been reconstructed that reproduces the phantom dark epoch with no 
need of any exotic fluid. In the same way, any given cosmology may be reconstructed. 



IV. UNIFIED INFLATION AND DARK ENERGY IN MODIFIED HORAVA-LIFSHITZ GRAVITY 

Let us consider here some viable F(R) gravities which admit the unification [l8| of inflation with late-time acceler- 
ation. In the convenient F(R) theory, a number of viable models (see [lj|, [2(|) which pass all local tests and are able 
to unify the inflationary and the current cosmic accelerated epochs have been proposed. Here we extend this class of 
models to the Hofava-Lifshitz gravity. We consider the action, 

F(R) = R + f(R), (35) 

where it is assumed that the term f(R) becomes important at cosmological scales, while for scales compared with the 
Solar system one the theory becomes linear on R. As an example, we consider the following function |19j 

- R n {aR n - 0) 
i{R) = l + 7 *» ' (36) 

where (a,/3,7) are constants and n > 1. This theory reproduces the inflationary and cosmic acceleration epochs in 
convenient F(R) gravity (see Ref. Q3), which is also the case in the present theory, as will be shown. During inflation, 
it is assumed that the curvature scalar tends to infinity. In this case the model (|35|. with (|36[) . behaves as 



Jim F(R) = aR n . (37) 

Then, by solving the FRW equation ([T5| , this kind of function yields a power-law solution of the type 

hi , , 2u(n - l)(2n - 1) 
Hit) = — , where h x = — ^ ^- '- . 38) 

This solution produces acceleration during the inflationary epoch if the parameters of the theory are properly defined. 
The acceleration parameter is given by - = hi(h\ — l)/£ 2 , thus, for hi > 1 the inflationary epoch is well reproduced 
by the model (|36[) . On the other hand, the function (|36[) has a minimum at Rq, given by 

Ro~(—) f , f(R)=0, /(£) = -2A~-£, (39) 



«7 j 7 



where we have imposed the condition /3 7 ja 3> 1. Then, at the current epoch the scalar curvature acquires a small 
value which can be fixed to coincide with the minimum (|3"9"|) . such that the FRW equations (fl"4"|) and (fT()|) yield 



(40) 
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which look very similar to the standard FRW equations in General Relativity, except for the parameter A. For first 



consideration of FRW eqs. in Hofava-Lifshitz gravity based on General Relativity see ref. (12|. As has been pointed 
out, at the current epoch the scalar R is small, so the theory is in the IR limit where the parameter A ~ 1, and the 
equations approach the usual ones for F(R) gravity. Hence, the FRW equations pt)|) reproduce the behavior of the 
well known ACDM model with no need to introduce a dark energy fluid to explain the current universe acceleration. 
As another example of the models described by (|35l) , we can considered the function (2^, [2l| , 



f{ R) = &-*>r +1 +%r +1 i Mil (41) 

fo + fi[(R-R ) 2n+1 + R 2 o n+1 \ h k + h [(R - Ro) 2n+1 + Rl n " 1 

This function could also serve for the unification of inflation and cosmic acceleration but, in this case, when one takes 
the limit R — > oo, one gets 

Jim F(R) = R - 2A, , where A; = 1/2 fx , (42) 

R— s-oc 

where the subscript i denotes that we are in the inflationary epoch. By inserting this into Eqs. (I14[) and (|16p . the 
FRW equations take the same form as in (|40|) . Then, for the function (|4Tj) the inflationary epoch is produced by an 
effective cosmological constant, which implies that the parameter A > 1/3, or the equations themselves will present 
inconsistencies, as it was discussed in the above section. For the current epoch, it is easy to see that the function (|4Tj) 
exhibits a minimum for R = Rq, which implies, as in the model above, an effective cosmological constant for late time 
that can produce the cosmic acceleration. The emergence of matter dominance before the dark energy epoch can be 
exhibited, in analogy with the case of the convenient theory. Hence, we have shown that the model (|41[) also unifies 
the cosmic expansion history, although with different properties during the inflationary epoch as compared with the 
model (J2U). This could be very important for the precise study of the evolution of the parameters of the theory. 

It is also interesting to explore the de Sitter solutions allowed by the theory iJJJ. By taking H(t) — Hq, the FRW 
equation (fTB|) . in absence of any kind of matter and with C — 0, reduces to 

= F(R ) - 6H 2 (1 - 3A + 3 f i)F'(R ) , (43) 

which reduces to an algebraic equation that, for an specific model, can be solved yielding the possible de Sitter points 
allowed by the theory. As an example, let us consider the model (|3l)f , where Eq. (|4"3"1) takes the form 

R n (ryR n - R) ®H 2 (-1 + 3A - 3/i) [l + na-fM 11 - 1 + i?"- 1 (2 7 i? - n0) + M* 1 ' 1 (~f 2 Rq + 2na) 
Ro + 01 °~ P) + l - = 1 = . (44) 



Here i?o = 3(1 — 3A + Q/i)Hq. By specifying the free parameters of the theory, one can solve Eq. ((44)) . which yields 
several de Sitter points, as the one studied above. They can be used to explain the coincidence problem, with the 
argument that the present will not be the only late-time accelerated epoch experienced by our Universe. In standard 
F(R), it was found for this same model that it contains at least two de Sitter points along the cosmic history (see 
22]). In the same way, the second model studied here (|41l) . provides several de Sitter points in the course of the 
cosmic history. Note that when fj, = A — |, Eq. (J43J) turns out to be much more simple, it reduces to F(Ro) = 0, 
where the de Sitter points are the roots. For example, for (l4"4l we have Ro(l + jRq) + Rq (o,Rq — j3) = 0, where the 
number of positive roots (de Sitter points) depends on the free parameters of the theory. 

Summing up, it has been here shown that, also in F(R) Hofava-Lifshitz gravity, the so-called viable models, as (|36|) 
or (l4Tj) , can in fact reproduce the whole cosmological history of the universe, with no need to involve any extra fields 
or a cosmological constant. We do not, of course, believe the complicated model (|36p or (|4"Tj) could be a final model. 
These models, at best, could be low energy effective theories coming from a more beautiful final theory. Even in the 
usual quantum field theories, the corresponding low energy effective theories are rather complicated due to the terms 
induced by the quantum effects. However, our main motivation to consider precisely them is related with the fact 
that they pass the known local tests as usual F(R) theories. 



V. NEWTON LAW CORRECTIONS IN F(R) GRAVITY 



As is well-known, modified gravity may lead to violations of local tests. We explore in this section how to avoid 
these violations of Newton's law. It is known that F(R) theories include scalar particle. This scalar field could give 
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rise to a fifth force and to variations of the Newton law, which can be avoided by a kind of the so-called chameleon 
mechanism [23[ . Note that scalars with time-dependent mass were also considered in 24] . 

In the original Hofava gravity, the projectability condition may cause problems with the Newton law [3l[ but the 
model without the projectability condition could be inconsistent for the F(i?)-model [ll[. The Newton law in the 



Hofava gravity may be restored by the "healthy" extension [3l| . In this section, we do not discuss the gravity sector 
corresponding to the Hofava gravity but we show that the scalar mode, which also appears in the usual F(R) gravity, 
can decouple from gravity and matter, and then the scalar mode does not give a measurable correction to Newton's 
law. 

To show this, we consider a function of the type (I55|) . and rewrite action (J7J as 

dtd 3 xy r ^)N\(l + f'(A))(R- A) + A + f(A)} , (45) 



where A is an auxiliary scalar field. It is easy to see that variation of action (|45|) over A gives A = R. Performing the 

glf = e-'t'~glf,wm<f> = l 



conformal transformation gff = e , with 0=1 ln(l + f'(A)), action (|4"5l) yields 



dtdPxy g( 



(46) 

where Kij is the extrinsic curvature given by as in <j3j> , K = g^^Kij, L(g( 3 \<fi) is the conformally transformed 
term in (|HJ|, and the scalar potential is 

A((j>)f'(A((j))) f{A{4>)) 

vw = iTTW ' (47) 

Note that, differently from the convenient F(R), in action (|46|) there is a coupling term between the scalar field (f> and 
the spatial metric <7« , which can thus be dropped, by imposing the following condition on the parameters 

H = X - | . (48) 

This condition also renders the theory power-counting renormalizable, for the same z as in the original Hofava model. 

Let us now investigate the term ©. As already pointed out, for local scales, where the scalar curvature is assumed 
to be very small, the theory enters the IR limit, where such term could be written as a spatial curvature, 

LW{g( 3 \<t>)~RW. (49) 

Then, corrections to the Newton law will come from the coupling that now appears between the scalar field and 
matter, which makes a test particle to deviate from its geodesic path, unless the mass of the scalar field is large 
enough (since then the effect could be very small). The precise value can be calculated from 

2 _l d?V(cp) _ l + r{A) A + f(A) 
m * "2 #2 - r{A) l + r{A) ■ 

In view of that, we can now analyze the models studied in the last section. We are interested to see the behavior 
at local scales, as on Earth, where the scalar curvature is around A = R ~ 10 _50 eV 2 , or in the solar system, where 
A = R ~ 10 _61 eV 2 . The function (|36| and its derivatives can be approximated around these points as 



f(R)^ r ^-^R^\ (51) 

7 7 7 



Then, the scalar mass for the model (|36l) is given approximately by the expression 

m l ~ 1^ u ' (52) 

^ n(n — l)a 

which becomes to 2 ~ io 50n - 100 e V 2 on Earth and m| - io 61 «- 12 2 j n tne Solar System. We thus see that, for n > 2, 
the scalar mass would be sufficiently large in order to avoid corrections to the Newton law. Even for the limiting case 
n = 2, the parameters j/a can be chosen to be large enough so that any violation of the local tests is avoided. 
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For the model (HU, the situation is quite similar. For simplicity, we impose the following condition 



fo < fi 



{R-Ro) 



2n+l 



n2n+l 

K 



fiR 



2n+l 



Then, function (|41[) and its derivatives can be written, for small values of the curvature, as 



f(R) 



fo 



fl f?R 2n+1, ' #R2(n+X) 

Using now the expression for the scalar mass (|50l) , this yields 

f?R 2n+3 



f(R) „ ( 2 "+ 1 )/0 ; r {R) „ _ 2 (2« + V( 2n + 2)/o 



f2 R 2n+3 



1 



2(2n + l)(n + l)/ /i 



(53) 



(54) 



(55) 



Hence, as R is very small at solar or Earth scales (~ 10 



-50 61 



eV 2 ), and Aj 



i 

2/i 



10 



20~38 



eV is the effective 



cosmological constant at inflation, which is much larger than the other term, it turns out that the second term on the 
rhs of ([55]) will dominate, and the scalar mass will take a value such as mi ~ 4- ~ 

10 20-38 eV 2 ; which ig large 

enougn 

as compared with the scalar curvature. As a consequence there is no observable correction to Newton's law. 

We have thus shown, in all detail, that the viable models (|36| and (|4Tj) do not introduce any observable correction 
to the Newton law at small scales. This strongly supports the choice of F(R) gravity as a realistic candidate for the 
unified description of the cosmological history. 



VI. FINITE-TIME FUTURE SINGULARITIES IN F(R) GRAVITY 

It is a well-known that a good number of effective phantom/quintessence- like dark energy models end their evolution 
at a finite-time future singularity. In the current section, we study the possible future evolution of the viable F(R) 
Hofava-Lifshitz gravity considered above. It has been already proven [5j that power-law F(R) Hofava-Lifshitz gravities 
may lead, in its evolution, to a finite-time singularity. In order to properly define the type of future singularities, let 
us rewrite the FRW Eqs. (fT4|) and (fT6|) in the following way 



3H 2 



3A- 1 



PcS ■ 



-3H 2 - 2H 



3A- 1 



PeS , 



(56) 



where 



1 



PeS 



PcS 



k 2 F'{R) 
1 



k 2 F'(R) 



—F(R) + 3(1 - 3A + 9n)H 2 F'(R) - QfiHF'(R) - 6/iHRF"(R) + k 2 p .„ 
F(R) - (6fiH + 3(1 - 3A + 9fj,)H 2 )F\R) - 2(1 - 3X)RF"(R) 



-2(i 



R 2 F {3) (R) + RF"(R) 



K 2 p r , 



(57) 



Then, using the expressions for the effective energy and pressure densities just defined, the list of future singularities, 
as classified in Ref. [25j], can be extended to F(R) gravity as follows: 

• Type I ("Big Rip"): For t — >• t s , a — > oo and p e s — > oo, \p\ — > oo. 

• Type II ("Sudden"): For t — > t s , a — ► a s and p c g — > p s , \p e s\ °°- 

• Type III: For t — > t s , a — > a s and p c g —> oo, \p e s\ ~> oo. 

• Type IV: For t —¥ t s , a —¥ a s and p c g — > p s , p c g — > p s but higher derivatives of Hubble parameter diverge. 

To illustrate the possibility of future singularities in viable F(R) gravity, we explore the model (I3"6")h which has 
been studied in Ref. [26[ for the convenient F(R) case, where it was shown that this model is non-singular, in the 
particular case n = 2. Also, it is known that for most models of F(R) gravity, the future singularity can be cured by 
adding a term proportional to R 2 (see Ref. [13, [28j]) or a non-singular modified gravity action (26|. However, this is 
not possible to do in the Hofava-Lifshitz gravity where, even for the simple model (|36[) with n = 2, a singularity can 
occur unless some restrictive conditions on the parameters are imposed. 



10 



In order to study the possible singularities that may occur from the above list, we consider a Hubble parameter 
close to one of such singularities, given by the following expression 

where ho and q are constant. Depending on the value of q, the Hubble parameter (|58p gives rise to a particular type 
of singularity. Thus, for q > 1, it gives a Big Rip singularity, for — 1 < q < a Sudden Singularity, for < q < 1 a 
Type III singularity and for q < —1, it will produce a Type IV singularity. Using the expression (|13[) with (|58[) , the 
scalar curvature can be approximated, depending on the value of q, as 

(t.-t)£ for « > 1 

3(l-3A+6M)fe|+6wfro f or g = 1 . (59) 
(t s -t)«+i t0r < 1 

The model (155)) which, as has been shown, unifies the inflationary and the dark energy epochs, makes the scalar 
curvature grow with time so that, close to a possible future singularity, the model can be approximated as F(R) ~ R n , 
where n > 1. For the case q > 1, it is possible to show, from the FRW Eqs. ([UP) and (fTlj)) . that the solution ((55)) is 
allowed for this model just in some special cases: (i) For n = 3/2 and /x = 2A — |, which contradicts the decoupling 
condition and the Newtonian corrections (where it was found that n > 2). (ii) When /i = and A = 1/3, which 
holds (|48|) and fixes the values of the parameters, in contradiction with the fact that fluctuations are allowed for them. 
For q — 1, the Hubble parameter (|58[) is a natural solution for this model, yielding 

ho , , 2u(n - l)(2n - 1) 
£T(i) = — — with ho = —, ~, r • 60) 

w t s -t -l + 3A-6/x + 2n(l-3A + 3^) V ' 

This model allows for the possibility of occurrence of a Big Rip singularity [29j , unless the parameters of the theory 
are fixed. As pointed out in Ref. [26j, for the case n = 2 the Big Rip singularity can be avoided in standard F(R) 
gravity, which can be easily seen by choosing A = \i = 1 in the solution (f6"0"]) . Nevertheless, in F(R) gravity, in order 
to avoid such singularity, the power n of the model has to be fixed to the value 

"=2 + 2-6A + 6^ (61) 

which can be interpreted as another constraint on the parameters of the theory (although, when the decoupling 
condition (|48[) is satisfied, no constraints can be imposed on n). In addition, note that ho = —h\ in (|38p . where we 
imposed h\ > 1 with the aim to reproduce the inflationary epoch, so that ho would be negative and the solution (|6"0")) 
would correspond to the Big Bang singularity. Thus, no future doomsday will take place. 

For the last case, when q < 1, we find that the Hubble parameter ([55)) can be a consistent solution when Eq. P5|) 
is satisfied and q = (n(l — n) — l)/n(n — 2), although if we impose n > 2, as it was found above, q < —1, which 
implies a future singularity of the type IV. The exception here is when n = 2, that avoids the occurrence of any type 
of singularity when q < 1. Nevertheless, even in the case that, for any reason, some kind of future singularity would 
be allowed in the model Q36p. one must also take into account possible quantum gravity effects, which may probably 
become important when one is close to the singularity. They have been shown, in phantom models, to prevent quite 
naturally the occurrence of the future singularity [30l | . 

In summary, we have here proven that the theory (|36[) can actually be free of future singularities, and thus that it 
can make a good candidate for the unification of the cosmic history. Note that, from the corrections to the Newton 
law studied in the previous section, as well as from imposing avoidance of future singularities, we can fix some of the 
parameters of the theory. The other constant parameters that appear in Q36p — which expresses the algebraic relation 
between the powers of the scalar curvature — could be fixed by comparing the cosmic evolution with the observed 
data, quite in the same way as has been successfully done for standard F(R) gravity in Ref. [22j]. 



VII. DISCUSSION 



In summary, we have here investigated the FRW cosmology of a non-linear modified Hofava-Lifshitz F(R) gravity 
theory which has a viable convenient counterpart. As shown in [J|, for a special choice of the parameters, the FRW 
equations are just the same in both theories, and that the cosmic history of the first literally coincides with the one 
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for the viable F(R) gravity. For a more general version of the theory, the unified description of the early-time inflation 
and late-time acceleration is proven to be possible too; however, the details of the cosmological dynamics are here 
different. Moreover, corrections to Newton's law are negligible for an extensive region of the parameter space. We 
have demonstrated the emergence of possible finite-time future singularities, and their avoidance, by adding extra 
higher-derivative terms which turn out to be qualitatively different, as compare with conventional F(R) cosmology. 

Using the approach of Refs. 0, one can construct more complicated generalizations of modified gravity with 
anisotropic scaling properties. For instance, one can obtain non-local or modified Gauss-Bonnet Hofava-Lifshitz 
gravities. Technically, this is of course a more involved task, as compared with the F(R) case. The corresponding 
cosmology can in principle be studied in the same way as in the present work, with expected qualitatively similar 
results, owing to the fact that the convenient cosmological model has been already investigated, in those cases. It is 
foreseen that the study of such theories will help us also in the resolution of the some problems for the theories with 
broken Lorentz symmetry. For instance, it is quite possible that a natural scenario of dynamical Lorentz symmetry 
breaking, which would reduce gravity to the Hofava limit, may be found in this way, which would be certainly 
interesting. 
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